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Beginning from the Ashtekar formulation of canonical general relativity, we derive a physical 
Hamiltonian written in terms of (classical) loop gravity variables. This is done by gauge-Hxing 
the gravitational fields within a complex of three-dimensional cells such that curvature and torsion 
vanish within each cell. The resulting theory is holographic, with the bulk dynamics being captured 
completely by degrees of freedom living on cell boundaries. Quantization is readily obtainable by 
existing methods. 


INTRODUCTION 


CONVENTIONS AND NOTATION 


The strong equivalence principle is one of the building 
blocks of general relativity (GR). It implies that an ap¬ 
propriate diffeomorphism can make spacetime flat at any 
point, while it remains generally curved elsewhere. This 
classical notion of locality however, is thought to extend 
only up to the Planck scale where various approaches 
to quantum gravity [1-6] suggest that a discrete space- 
time structure emerges. The main challenge facing all of 
these programs is to develop a complete, self-consistent 
quantum dynamics that agrees with GR in the classical 
limit. This is the ‘bottom-up’ view. Here we take the 
‘top-down’ view, beginning from the continuous theory 
to first derive a discrete, dynamical theory of GR before 
constructing the quantum theory [7]. The main idea is 
to extend local flatness to small but finite ‘cells’ of space, 
over small (continuous) time intervals, and to impose this 
expicitly within the canonical formalism of GR. 

We define spatial hypersurfaces in terms of a count¬ 
able but arbitrary number of three-dimensional cells, con¬ 
straining the geometrical fields defined within each cell 
to be intrinsically and extrinsically flat. This method 
yields an invertible map from the continuous to a dis¬ 
crete phase space written in terms of holonomy-flux vari¬ 
ables for loop gravity [8]. These parameterize equivalence 
classes of continuous, piecewise-flat spatial hypersurfaces 
in terms of finite degrees of freedom (DOF), in a way 
which reproduces the full class of spatial hypersurfaces 
as the number of cells grows large [9]. On the reduced 
phase space, the constraints within each cell are identi¬ 
cally satisfied leaving a discrete theory written in terms 
of holonomies and fluxes on cell boundaries. Maintaining 
a continuous geometry dynamically completely fixes the 
gauge thereby providing a physical Hamiltonian. The 
result is a discrete theory of GR suitable for canonical 
quantization in terms of the operators and Hilbert space 
of loop quantum gravity (LQG). 


We use su(2) basis elements r® (for i = 1,2,3) which 
are given by —f/2 times the Pauli matrices. We define 
a trace notation such that Tr(rV-l) = —2trace(rV'l) = 
J*-!, and a bracket [•,•] which implies taking the su(2) 
commutator and wedge product between entries. Ele¬ 
ments V G su(2) are written in bold font where v = 
uV®, and su(2) indices are all written ‘up’ for conve¬ 
nience (since the Gartan-Killing metric is the Kronecker 
delta). Elements of su(2) represent vectors in BR^ via 
u* = Tr(r®t>), and we define the modulus |i;| = -\/Tr(i;t>). 


CONTINUOUS HAMILTONIAN 

A Hamiltonian for GR may be written in terms of 
the Ashtekar variables {A, E) [10]. The momentum two- 
form, or ‘electric field’, E describes the intrinsic geome¬ 
try of spatial hypersurfaces, being given in terms of the 
frame-field e according to: E = i[e, e]. The configura¬ 
tion variable is the connection one-form A. It contains 
both intrinsic and extrinsic information, being given by: 
A = r -|- iK, where K is the extrinsic curvature and F 
is the (torsion-free) Levi-Givita connection. 

At this point, the fields are defined on a simply con¬ 
nected three-manifold c, which has a closed boundary dc. 
The fields (A, E) together with the manifold c define a 
spatial hypersurface of spacetime. We begin by analyz¬ 
ing {A, E) within a single cell c, before discussing how 
to glue multiple cells together later on. 

The canonical action I = Ul — H is composed of the 
symplectic term Q and the Hamiltonian H. The sym- 
plectic term is given by: 

n = i J^Tr Aae, (1) 

which implies that { A, E) form a conjugate pair of phase 
space variables. The Hamiltonian is a sum of the scalar, 
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vector and Gauss constraints: 


^[Af] + G[|] = Tr^(7V + |)eAF, (2a) 
G[A] = Tr / XdAE, (2b) 


where N, ^ and A are Lagrange multipliers and F = 
dA+^[A,A] is the curvature of the Ashtekar connec¬ 
tion. The scalar constraint generates dynamics, the vec¬ 
tor constraint generates diffeomorphisms and the Gauss 
constraint generates SU(2) gauge transformations. The 
fields {A, E) and Hamiltonian (2) define GR with van¬ 
ishing cosmological constant. 

In order that the variational principal be well-defined, 
we add the following boundary terms to the Hamiltonian: 


Sb[N]+Vb[^] = Tt {N + ^)eAA, (3) 
J dc 

Gb[A] = -Tr [ XE. (4) 

J dc 

and place boundary conditions on the electric field as 
follows. Points in c are labeled by a chart x, and for sim¬ 
plicity we choose this to be a set of cartesian coordinates 
in which c takes the form of a tetrahedron. This shape 
has no meaning - the physically meaningful geometry is 
that described by E{x) which defines the spatial metric. 
We want this physical geometry to be a tetrahedron as 
well, but it is important to keep in mind that the tetra¬ 
hedron defined by E{x) is generally different than the 
one described by the chart x. In order that E describes 
a tetrahedron, we impose the boundary condition that E 
is constant along each face. The value of this constant is 
left free, which leaves three DOF in E per face and allows 
the faces to move. This boundary dynamics is our main 
result and is developed in the remainder of the article. 


BOUNDARY THEORY 

We now present a gauge reduction in which the canoni¬ 
cal action vanishes within c and leaves a boundary theory 
on dc. After studying a single three-manifold c, we sub¬ 
sequently consider a collection of these cells which are 
glued together to form a larger, piecewise smooth three- 
manifold that is not simply connected. In the end we 
arrive at a physical Hamiltonian describing the full dy¬ 
namics of such data, parameterized completely by DOF 
on the boundaries dc. 


Kinematics 

We impose conditions such that {A, E) describe an 
intrinsically and extrinsically flat geometry, as in the ‘flat 
cell gauge’ of [8] which was further studied in [11]. In 


terms of our phase space, this flatness implies [8] that the 
connection is flat F — 0 and its torsion vanishes: T = 
de + [A,e] = 0. With these conditions the fields satisfy 
the constraints (2a, 2b) and therefore provide initial data 
for GR. 

The solution to F = T = 0 is given in terms of 
an algebra-valued ‘coordinate’ function y and a group¬ 
valued ‘rotation’ function a: 


A = a ^da, e = a ^dya. (5) 

The rotation function a(x) is a holonomy defining parallel 
transport from some basepoint to any other point x G c 
along an arbitrary path [21]. The coordinate function y 
defines the intrinsic geometry of c. 

In order to develop an invertible map between the con¬ 
tinuous variables (A, E) and the discrete boundary vari¬ 
ables defined below, we choose a specific representative 
of the solutions (5) by defining (a,y) as follows. We 
choose each function y^{t,x) to be linear in x so that 
the coordinate function provides a different, time depen¬ 
dant chart in which c appears as a tetrahedron. We also 
fix y = 0 at the barycentre so that the tetrahedron is 
uniquely defined in terms of y-coordinates. Notice that 
c may change shape according to y, while it is static in 
terms of the chart x. 

To define the rotation function, we must first define 
some further structure. Gonsider a subdivision of c into 
four tetrahedron-shaped regions rf, where each face f is 
the base of a sub-tetrahedron, and the barycentre n is at 
its apex [22]. Now consider a one-parameter family of two 
surfaces labeled by constant values of (p{x) such that (p = 
0 is the union of ‘lateral’ faces (the faces which intersect 
at n) and = 1 is the face f. Notice however that <p{x) 
is multiply-defined on edges (or ‘bones’, in analogy with 
Regge geometry [12]) of the tetrahedron. 

In order to make p well-defined, we excise the bones 
as follows. Gonsider a cylinder of radius e with a bone 
b along the axis, parameterized by angular and longitu¬ 
dinal coordinates {4>,z). The boundary is defined in the 
limit of e —>■ 0 by keeping the angular points distinct [23]. 
This ‘smears’ out the values of ip over the boundary 9b so 
that it is now well-defined, varying as one travels around 
9b. 

The bones are to be seen as topological defects, and 
must remain one-dimensional in terms of y-coordinates. 
This requires the following boundary condition [14, 15]: 


^vy\db ~ 


( 6 ) 


Now we are ready to define the rotation function. Us¬ 
ing a bump function f{p) and a ‘twist’ parameter yf(t) 
we define a{x) for all x in rf as: 


a(t, p) = exp 



fip) 


( 7 ) 


for any point in a region rf. Notice path-ordering is not 
required since pf is constant in x. Due to the properties 
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of a bump function, one can see that the connection is 
smooth everywhere in c and vanishes on faces. Since the 
chart X is time-independent, the entire dynamics of a is 
contained within the twist parameters. 

The excision defined above is key for two reasons: 1) It 
allows curvature to be supported on bones while the fields 
{A, E) remain well-defined everywhere else; 2) treating 
these sources of curvature as topological defects prevents 
the curvature from spreading out as time evolves, so that 
the construction is preserved dynamically. 

Reduction 

The discrete boundary theory is obtained by using (5) 
as gauge conditions in the continuous theory. Maintain¬ 
ing these conditions throughout the evolution fixes the 
Lagrange multipliers, up to undetermined constants, in 
terms of (a, y) and (a,y). The fully-determined, non¬ 
constant parts describe the local diffeomorphisms and 
rotations necessary to preserve the form of (a, y) given 
above. The undetermined components are written as N, 
a~^\a and where an overbar denotes a constant; 

these determine global transformations over c and will be 
fixed in the subsequent section. 

Let us now impose (5) within the theory by direct sub¬ 
stitution in the canonical action. One finds that the 
terms (2a, 2b) each vanish identically, while the sym- 
plectic term vanishes in the bulk but leaves behind a 
boundary contribution [8, 9]: 

12 =-iTr^V^/ifXf. (8) 

f 

The four pairs {hf,Xf) are a discrete set of boundary 
variables: 

h^ = a~^{x), Xf = i^[dy,dy], (9) 

for any point x on f. Each holonomy /if defines parallel 
transport from a face f to the basepoint n. Each ‘flux’ Xf 
defines a vector normal to f with a modulus equal to the 
area. The geometry of the tetrahedron c is determined 
by these four fluxes. 

The discrete variables (/if,Xf) parameterize a phase 
space [13] on each face with the following Poisson brack¬ 
ets: 

{/if, /le} = 0, {x/, X/;} = i,5fpeb'=X/=, (10) 

{/if,X/,} =/(5ff/r*/if, (11) 

where f and f' are any two faces on the boundary dc. 

Let us now turn to the boundary terms. The fully- 
determined parts of the multipliers lead to terms which 
do not contribute to the discrete theory. By direct cal¬ 
culation, the undertermined components yield for (4): 

GB[A] = -TrA^Xf. (12) 

f 


This is the well-known discrete Gauss (or closure) con¬ 
straint which tells us that the fluxes add up to zero, as 
required for a closed two-surface. 

Looking at the boundary terms (3), we find that these 
vanish on faces because A is constant along each face. 
On the two-surface i9b surrounding each bone however, 
we find: 

SB[N] + VB[i]=Ty{N + Oj2 E (13) 

f bipbnf/o 

The first sum is over faces, and the second is over bones 
which intersect a face. This simple form is obtained using 
that: 

/ dip d^aa~^ = pf, (14) 

J ^brirf 

and defining a vector: 

Qb ■■= [ dz dzV, (15) 

J 5bnrf 

which gives the length and orientation of the bone. 

In order to make use of (13) in the boundary theory, 
we need to write it in terms of the phase space variables 
(/if,Xf). This is possible because of the way in which 
the fields (a,y) have been defined. The Qb define edges 
of a tetrahedron, and these can be written explicitly in 
terms of fluxes Xf using simple geometrical arguments. 
For the twists, note that an element of the group SU(2) 
is the exponential of an su(2) algebra element /if = e^', 
where pf defines an angle jpfj and axis of rotation. This 
map is invertible for 0 < jpfj < 27r. 

The canonical theory for {A, E) in c has been ex¬ 
pressed entirely in terms of four sets of boundary vari¬ 
ables (/if,Xf). This discrete theory has the following 
Hamiltonian: 

iLc = Sb [X] + Vb ® + Gb [A], (16) 

where each term is given above in (12, 13); the vectors 
Qb are explicit functions of fluxes, and the twists pf are 
explicit functions of the holonomies. At any given time, 
knowing the holonomies and fluxes allows one to unam¬ 
biguously determine helds {A,E) from (5), using that 
y* are cartesian coordinates for c with y* = 0 at the 
barycentre, and a is given by (7). Therefore, evolution in 
the bulk is completely determined by (16) in this reduced 
theory. 

For a single manifold c, the holographic correspon¬ 
dence described above holds for any choice of multipliers 
(X, 4, A) which are spatially constant in c. In the next 
section we glue many c’s together to form a larger space. 
Requiring that the resulting geometry is continuous at 
all times fixes the value of each multiplier, within each c. 

In the boundary theory, variables (/if,Xf) are associ¬ 
ated to each of the four faces in dc giving 6 x 4 = 24 
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(complex-valued) DOF. These are subject to the (first 
class) closure constraint which removes six DOF and 
leaves 18 phase space DOF to describe the evolution. 
It is interesting that this is the same number of DOF 
that the Ashtekar variables possess at each point in the 
continuous theory, before taking gauge freedom into ac¬ 
count. This suggests an analogy between points in the 
continuous theory and cells in the discrete theory. 

GLUING 

We now generalize to the case where c is but one of 
many ‘cells’ in a CW complex [16], each constructed in 
the manner described above. In order that the fields 
{A, E) describe a continuous geometry throughout the 
complex, certain gluing conditions must hold so that the 
triangles shared by adjacent tetrahedra match up. These 
are gauge conditions, and as we will see, maintaining 
these conditions dynamically fixes a unique solution for 
{N, A). We refer the reader to [15] for a more detailed 

description (with illustrations) of a similar gluing in two 
dimensions, with the caveat that there are important dif¬ 
ferences due to the trivial nature of three-dimensional 
gravity. 

To keep the presentation simple, let us consider a CW 
complex with only two cells c and c', where each face f S 
dc is identified with a face f' S dc'. The results presented 
below generalize to an arbitrary number of cells. 

Within each cell the fields satisfy F = T = 0, and 
the solution is given in terms of local functions (ac,yc) 
and (ocsyc')- In order that the geometry is continuous, 
both cells must agree on the value of E at each shared 
face. From the form of E given by (5), this implies that 
differentials of the coordinate functions must agree, up 
to a rotation, as one approaches each face from either 
side. This provides four conditions, one for each face. 
Note that A goes to zero smoothly as one approaches a 
face from either side, so that no gluing is required for a 
continuous connection. 

These conditions imply that the flux Xf as seen from 
c is related to the flux Xfi seen from c' by: 

Xf = -gfXf,g-\ (17) 

This ‘gluing’ condition fixes the areas of f and f' to be 
equal, i.e. |Xf| = |Xf/|. The relative orientation is a free 
variable gt = describing parallel transport from c 

to c' through the face f = f'. 

Although the areas of each identified pair (f, f') are 
now set equal, the shapes are generally different, so at 
this point we have a (discontinuous) ‘twisted’ geometry 
[17]. To obtain a continuous geometry we must also im¬ 
pose shape-matching conditions [18] which identify an¬ 
gles shared by f and f' so that the triangles represented 
by each face are the same. With these conditions and 
(17), the fluxes now define a spatial Regge geometry. 


Let us count the number of shape matching conditions 
needed on this two-cell CW complex. A single tetrahe¬ 
dron is defined by six parameters, so six conditions in to¬ 
tal are needed to give each tetrahedron the same shape. 
We have given four above in (17) to match the areas, 
so we require two shape-matching conditions to obtain a 
continuous geometry across each face in the two-cell CW 
complex. 

With the gluing and shape-matching conditions we 
have obtained a continuous geometry everywhere on the 
CW complex. Each tetrahedron-shaped cell is parame¬ 
terized by {hf, Xf) subject to these conditions and the 
closure constraint. The dynamics is generated by a 
Hamiltonian (16) local to each cell, which can be written 
explicitly in terms of the holonomy-flux variables defining 
each cell. 

The number of DOF remaining after gluing are counted 
as follows. Consider two (so far unglued) cells for a total 
of eight faces, with a six parameter phase space {hf,Xf) 
on each face. We must then subtract: two closure con¬ 
straints (one for each cell); four gluing conditions (one for 
each pair of faces); two shape-matching conditions (one 
for each cell). The closure constraints and gluing condi¬ 
tions each contain three equations while shape-matching 
is a single equation; each equation removes a pair of phase 
space DOF. Considering all of this, the total DOF are 
6x (8 — 2 — 4) — 2x2 = 8. For an arbitrary number 
of cells one finds four DOF per cell, the same number of 
physical DOF as GR has per point. 

The previous paragraph implies that gluing cells to¬ 
gether to form a continuous geometry completely reduces 
the theory. The covariantly-constant parts (fVc,^c,Ac) 
of the Lagrange multipliers must then have been fixed 
by the gluing and shape-matching conditions. The equa¬ 
tions which do so come from the time derivatives of these 
conditions, using the Hamiltonian equations of motion 
for Xf and hf. In the two-cell complex these total 14 
equations, linear in {Nc,^c, ^c), which exactly fix these 
14 multipliers. For an arbitrary number of cells, one ob¬ 
tains seven equations per cell to fix seven multipliers per 
cell. Note that these solutions are highly non-local, in¬ 
volving data from every cell in the complex. Notice also 
that after gluing, the analogy still holds between points 
in the continuous theory and cells in the discrete theory. 

DISCUSSION 

We have obtained a discrete Hamiltonian for gravity 
from a reduction of the continuous theory. Beginning 
from the Ashtekar variables {A, E) for canonical GR, we 
reduced to fields which are flat and torsion free in a piece- 
wise manner, satisfying F = T = 0 within each cell of an 
GW complex. For such a geometry, curvature and torsion 
have support only upon the one-skeleton, i.e. the bones. 
Gluing the cells together to form a continuous spatial 
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geometry, and maintaining continuity at all times, gives 
conditions which completely fix the gauge. This yields 
a physical Hamiltonian (16) which generates dynamics 
in the discrete phase space. The resulting theory is not 
an approximation, but rather an exact description of dy¬ 
namical, piecewise-flat and piecewise-torsion-free geome¬ 
tries. 

The discrete theory is written within the holonomy- 
fiux phase space of loop gravity, from which the Hilbert 
space and operators of LQG [3] can be constructed. 
These well-established techniques can be readily applied 
here — the new feature is a Hamiltonian operator com¬ 
ing from the quantization of (16). To help develop this 
operator, the discrete theory serves as a valuable consi- 
tency check since the quantum dynamics must reproduce 
the classical dynamics in the h ^ 0 limit. Considering 
then that the discrete phase space in the continuum limit 
(when the number of cells grows large) describes the usual 
spatial geometries for canonical GR [9], quantization of 
the discrete theory is a promising avenue toward quan¬ 
tum gravity. 
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